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MARK WATKINS 

Abstract. We give some heuristics for counting elliptic curves with certain 
properties. In particular, we re-derive the Brumer-McGuinness heuristic for 
the number of curves with positive/negative discriminant up to X, which is 
an application of lattice-point counting. We then introduce heuristics (with 
refinements from random matrix theory) that allow us to predict how often 
we expect an elliptic curve E with even parity to have L{E, 1) = 0. We find 
that we expect there to be about ciX^^/^'*(log X)'^/* curves with |A| < X 
with even parity and positive (analytic) rank; since Brumer and McGuinness 
predict cX^/^ total curves, this implies that asymptotically almost all even 
parity curves have rank 0. We then derive similar estimates for ordering by 
conductor, and conclude by giving various data regarding our heuristics and 
related questions. 

1. Introduction 

We give some heuristics for counting elliptic curves with certain properties. In 
particular, we re-derive the Brumer-McGuinness heuristic for the number of curves 
with positive/negative discriminant up to X, which is an application of lattice-point 
counting. We then introduce heuristics (with refinements from random matrix 
theory) that allow us to predict how often we expect an elliptic curve E with even 
parity to have L{E, 1) = 0. It turns out that we roughly expect that a curve with 
even parity has L{E, 1) = with probability proportional to the square root of its 
real period, and, since the real period is very roughly this leads us to the 

prediction that almost all curves with even parity should have L{E, 1) 7^ 0. By 
the conjecture of Birch and Swinncrton-Dycr, this says that almost all such curves 
have rank 0. 

We then make similar heuristics when enumerating by conductor. The first task 
here is simply to count curves with conductor up to X, and for this we use heuristics 
involving how often large powers of primes divide the discriminant. Upon making 
this estimate, we are then able to imitate the argument we made previously, and 
thus derive an asymptotic for the number of curves with even parity and L{E, 1) = 
under the ordering by conductor. We again get the heuristic that almost all curves 
with even parity should have L{E, 1) 7^ 0. 

We then make a few remarks regarding how often curves should have nontriv- 
ial isogenics and/or torsion under different orderings, and then present some data 
regarding average ranks. We conclude by giving data for Mordell-Weil lattice dis- 
tribution for rank 2 curves, and speculating about symmetric power L-functions. 



2. The Brumer-McGuinness Heuristic 

First we re-derive the Brumer-McGuinness heuristic 3 for the number of elliptic 
curves whose absolute discriminant is less than a given bound X] the technique here 
is essentially lattice-point counting, and we derive our estimates via the assumption 
that these counts are well-approximated by the areas of the given regions. 
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Conjecture 2.1. [Brumer-McGuinness] The number A±(X) of rational elliptic 
curves with a global minimal model ( including at oo) and positive or negative dis- 
criminant whose absolute value is less than X is asymptotically A±{X) ~ ^^y-'^^^^^ 

where a± = ^J^^^^. 

As indicated by Brumcr and McGuinncss, the identity a- = ^/3a+ was already 
known to Legendre, and is related to complex multiplication. These constants can 
be expressed in terms of Beta integrals B{u, v) = /g x^~^{l — x)""! dx = ^-^^u+v) 
as a+ = iB(l/2, 1/6) and a. = S(l/2, 1/3). 

Recall that every rational elliptic curve has a unique integral minimal model 
+ aixy + a^y = x^ + a2X^ + 040; + 05 with ai, 03 G {0, 1} and \a2\ < 1. Fix one 
of the 12 choices of (01,02,03). Since these are all bounded by 1 the discriminant 
is thus approximately — 64a4 — 432a|. So wc wish to count the number of (0,4, ar,)- 
lattice-points with |64a4 + 432ag| < X, noting that Brumer and McGuinness divide 
the curves according to the sign of the discriminant. The lattice-point count for 
Oi = 02 = 03 = is given by 



1 + 



We estimate this lattice-point count by the integral j J^j du^duQ for the region U 
given by |64u4 + 432ug| < X. After splitting into two parts based upon the sign of 
the discriminant and performing the W4-integration, we get 



(64)1/3 



f 

Jo 



(-432^2)1/3 _ (_x - 4321*2)1/3 



dU6 + 



(64)1/3 



[X - 432«2)V3 _ (_432u2)i/3 



due, 



where the factor of 2 comes from the sign of uq. Changing variables uq = w^Xj^ZI 
and multiplying by 12 for the number of (oi, 02, a3)-choices we get 



^5/6 



\v? + 1)1/3 




(«;2)i/3 - iw" - 1) 



1/3 



dw. 



These integrals are probably known, but I am unable to find a reference. The 
first integral simplifies^ to 1 ^^^^ 1 ^ Z^i^/'^^ 1/6), while the second becomes 



I J°°^ ^^+1 ~ 1-^(1/-^' 1/3)- This counts all models of curves; if wc eliminate non- 
minimal models, for which we have p^|c4 and p^^l^ fgp some prime p, we expect to 
accrue an extra factor^ of C(10)- Ftom this, we get the conjecture of Brumer and 
McGuinness as stated above. 



IAs N. D. Elkics indicated to us, wc can write /(a) = [(t^ -|-a)^/^ - (t^)!/-^] dt, differentiate 
under the integral sign, then substitute + a = ax'^, and finally re-integrate to obtain /(I). 

^Note that some choices of (ai, 02, as) necessarily have odd discriminant, but the other choices 
compensate to give the proper Euler factors at 2 (and 3). 
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3. Counting curves of even parity whose central X-value vanishes. 

See |2HI §15-16] for definitions of tlie conductor N and L-function L{E,s) of an 
elliptic curve E. Since rational elliptic curves are modular, we have that the com- 
pleted L-function A(£', s) = r(s)(%/]V/27r)"i(£', s) satisfies A(£', s) =±A(£',2-s). 
When the plus sign occurs, we say that E has even parity. 

We now try to count elliptic curves E with even parity for which L{E, 1) = 0. 
Throughout this section, E shall be a curve with even parity, and we shall order 
curves by discriminant. Via the conjectural Parity Principle, we expect that, under 
any reasonable ordering, half of the elliptic curves should have even parity; in 
particular, we predict that there are asymptotically A±{X)/2 curves with even 
parity and positive/negative discriminant up to X. 

Our main tool shall be random matrix theory, which gives a heuristic for pre- 
dicting how often L{E, 1) is small. We could alternatively derive a cruder heuristic 
by assuming the the order of the Shafarevich-Tate group is a random square integer 
in a given interval, but random matrix theory has the advantage of being able to 
predict a more explicit asymptotic. Our principal heuristic is the following: 

Heuristic 3.1. The number R{X) of rational elliptic curves E with even parity 
and L{E, 1) = and absolute discriminant less than X is given asymptotically by 
R{X) ^ cX^^/^'^(logX)^/^ for some computable constant c > 0. 

In particular, note that we get the prediction that almost all curves with even 
parity have L(E, 1)^0 under this ordering. 

3.1. Random matrix theory. Originally developed in mathematical statistics by 
Wishart [S^ in the 1920s and then in mathematical physics (especially the spectra of 
highly excited nuclei) by Wigner j3Sl, Dyson, Mehta, and others (particularly [2]), 
random matrix theory |23j has now found some applications in number theory, the 
earliest being the oft-told story of Dyson's remark to Montgomery regarding the 
pair-correlation of zeros of the Riemann C- function. Based on substantial numerical 
evidence, random matrix theory appears to give reasonable models for the distri- 
bution of L-values in families, though the issue of what constitutes a proper family 
is a delicate one (see particularly §3] , where the notion of family comes from the 
ability to calculate moments of L- functions rather than from algebraic geometry). 

The family of quadratic twists of a given elliptic curve E : y'^ = + Ax -I- i? is 
given by Ed : y^ = x^ + Ad^x + Bd'^ for squarefree d. The work (most significantly a 
monodromy computation) of Katz and Sarnak jl7| regarding families of curves over 
function fields implies that when we restrict to quadratic twists with even parity, 
we should expect that the L-functions are modelled by random matrices with even 
orthogonal symmetry. Though we have no function field analogue in our case, we 
brazenly assume (largely from looking at the sign in the functional equation) that 
the symmetry type is again orthogonal with even parity. What this means is that 
we want to model properties of the L-function via random matrices taken from 
S0(2M) with respect to Haar measure. Here we wish the mean density of zeros of 
the L-functions to match the mean density of eigenvalues of our matrices, and so, as 
in US], we should take 2M « 21ogA^. We suspect that the L-value distribution is 
approximately given by the distribution of the evaluations at 1 of the characteristic 
polynomials of our random matrices. In the large, this distribution is determined 
entirely by the symmetry type, while finer considerations are distinguished via 
arithmetic considerations. 
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With this assumption, via the moment conjectures of ^H] and then using MeUin 
inversion, as t ^ we have (see (21) of [J]) that 

(1) Pmh[L{E, 1) < t] - a{E)t^l'^M^/^. 

This heuristic is stated for fixed M « logiV, but we shall also allow M oo. 
It is not easy to understand this probability, as both the constant a{E) and the 
matrix-size M depend on E. We can take curves with e*^ < N < e^^'^^ to mollify 
the impact of the conductor, but in order to average over a set of curves, we need to 
understand how a{E) varies. One idea is that a{E) separates into two parts, one 
of which depends on local structure (Frobenius traces) of the curve, and the other 
of which depends only upon the size of the conductor TV. Letting G be the Barnes 
G-function (such that G{z + 1) = r{z)G{z) with G(l) = 1) and M = [logTVj we 
have that 

a{E) ^ aniM) ■ aAiE) 

with afl(M) ^Ar^ 2i/**G(l/2)7r-i/4 as M oo and 
(2) 

^ n-<.. ^ n(- in^) {y-^^ 

where Lp{X) = {l — UpX+pX'^)^-^ whcnpj A and Lp{X) = {1 — apX)^^ otherwise; 
see (10) of evaluated at A: = —1/2, though that equation is wrong at primes that 
divide the discriminant — see (20) of |H], where Q should be taken to be 1. Note 
that the Sato- Tate conjecture 1311 implies that is p on average, and this implies 
that the above Euler product converges. 

3.2. Discretisation of the L-value distribution. For precise definitions of the 
Tamagawa numbers, torsion group, periods, and Shafarevich-Tate group, see j28| . 
though below we give a brief description of some of these. We let Tp{E) be the Tam- 
agawa number of E at the (possibly infinite) prime p, and write t(E) = Y[p TpiE) 
for the Tamagawa product and T{E) for the size of the torsion group. We also 
write fti-c{E) for the real period and nian(-E) for the size of the Shafarevich-Tate 
group when L{E, 1) 7^ 0, with m^niE) = when L{E, 1) = 0. 

We wish to assert that sufficiently small values of L{E, 1) actually correspond 
to L{E, 1) — 0. We do this via the conjectural formula of Birch and Swinnerton- 
Dyer (Ij, which asserts that 

L{E,i) = n,.,{E)-^^-m,4E). 

Our discretisation'^ will be that 

L{E,l)<n,,iE)-^^ implies L{E,l)^0. 

Note that we are only using that lUan takes on integral values, and do not use the 
(conjectural) fact that it is square. 



■^The precision of this discretisation might be the most-debatable methodology we use. Indeed, 
we are essentially taking a "sharp cutofT" , while it might be better to have a more smooth transition 
function. For this reason, we do not specify the leading constant in our final heuristic. 
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Using we estimate the number of curves with positive (for simphcity) dis- 
criminant less than X and even parity and L{E^ 1) = via the lattice-point sum 



W{X)= a«(MW(i?)yM|^.M3/B. 

0<c|-c2<172SX 

We need to introduce congruence conditions on C4 and cg to make sure that they 
correspond to a minimal model of an elliptic curve. The paper ^01 uses the 
work of Connell in a different context to get that there are 288 classes of 
(c4 mod 576, cg mod 1728) that can give minimal models, and so we get a factor 
of 288/(576- 1728), assuming that each congruence class has the same impact on all 
the entities in the sum. Indeed, this independence (on average) of various quantities 
with respect to C4 and cg is critical in our estimation of W{X). There is also the 
question of non-minimal models,^ from which we get a factor of 1/^(10). 

Guess 3.2. The lattice-point sum W{X) can he approximated as X 00 by 

= (571^728) CM ■ "^'^"^^ • / / ■ (l^S^)^^^-^"^^"- 

1< T7^X^ <^ 



Here d/j is the limit I'^I^GiXjl^'K^^I^ of aji{M) as M — > 00, while is a suitable 
average of the arithmetic factors aA{E), and (3^ is the average of the square root 
of the Tamagawa product. We have also approximated log N « log A and assumed 
the torsion is trivial; below we will give these heuristic justification (on average). 
Note that everything left in the integral is a smooth function of U4 and uq . 

We shall first evaluate the integral in W(X) given these suppositions, and then 
try to justify the various assumptions that are inherent in this guess. ^ 

3.3. Evaluation of the integral. Write E as — Ax^ — (c4/12)a; — C6/216, and 
put ei > 62 > 63 for the roots of the cubic polynomial on the right side. We have 

l/rirc = agm(Vei - 62, \/ei - 63) /tt. 



We also have that (ei — e2)(ei — e3)(e2 — 63) — ^A/16 from the formula for the 
discriminant. We next write ei — 62 = A^/^A and 62 — 63 = A^/^/i so that we have 
^A(A-(-Ai) = 1/4, while d = ^{n+2\), e^ = ^{^-X), and 63 = -Al^(2^+A). 
Thus we get 

Ai/2 

-C6/864 = -616263 = ^(m + 2A)(/i - A)(2/i + A) 

and 

Ai/3 

-C4/48 = 6162 + 6163 + 6263 = ^(/^^ + A^ + A^). 

Changing variables in the VF-integral gives a Jacobian of 432/ A^^^ y/JJ + ^ so that 

(logA)3/8 dfldA 



WiX) 



^Ai/12 agm(\/A, VA+TI) ^^/Vm^ + m' 



^At p = 2, 3, non-minimality occurs when C4/p'^ and c^/p^ satisfy the congruences 
^Note that our methods do not readily generalise to higher rank, as there is no apparent way 
to model the heights of points (and thus the regulator). 
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where A — ( + M ^ Thus the variables are nicely separated, and since 

the /z-integral converges, we do indeed get W{X) - cX^^/'^^{\ogXf/^. A simi- 
lar argument can be given for curves with negative discriminant. This concludes 
our derivation of Heuristic 13.11 and now we turn to giving some reasons for our 
expectation that the arithmetic factors can be mollified by taking their averages. 

3.4. Expectations for arithmetic factors on average. In the next section we 
shall explain (among other things) why we expect that log N « log A for almost all 
curves, and in sectional we shall recall the classical parametrisations of Xi{N) due 
to Fricke to indicate why we expect the torsion size is 1 on average. Here we show 
how to compute the various averages (with respect to ordering by discriminant) of 
the square root of the Tamagawa product and the arithmetic factors aA{E). 

For both heuristics, we shall make the assumption that curves satisfying the 
discriminant bound |A| < X behave essentially the same as those that satisfy 
|c4| < X^^^ and \ce\ < X^l"^ . That is, we approximate our region by a big box. We 
write D for the absolute value of A. First we consider the Tamagawa product. 

We wish to know how often a prime divides the discriminant to a high power. 
Fix a prime p> 5 with p a lot smaller than X^^^. We can estimate the probability 
that I A by considering all p'^'^ choices of C4 and cg modulo p'', that is, by counting 
the number of solutions S{p'') to c\ — c\ — 1728 A = (mod p^). This auxiliary 
curve C4 = Cg is singular at (0, 0) over Fp, and has (p— 1) non-singular Fp-solutions 
which lift to p^~^{p — 1) points modulo p^ . 

For p^ sufficiently small, our (04, C6)-region is so large that we can show that the 
probability that p'^jA is S{p^)/p'^''. We assume that big primes act (on average) in 
the same manner, while a similar heuristic can be given for p = 2, 3. Curves with 
p^|c4 andp^jcg will not be given by their minimal model; indeed, we want to exclude 
these curves, and thus will multiply our probabilities by Kp = (1 — 1 /_p^°)~^ to make 
them conditional on this criterion. For instance, the above counting of points says 
that there is a probability of {p^ — p)/p'^ that p \ D, and so upon conditioning upon 
minimal models we get Kp(l — 1/p) for this probability. 

What is the probability Pm{p,k) that a curve given by a minimal model has 
multiplicative reduction at p > 5 and p'^WD for some fc > 0? In terms of Kodaira 
symbols, this is the case of Ifc. For multiplicative reduction we need that p \ 04, cq. 
These events are independent and each has a probability (1 — 1/p) of occurring. 
Upon assuming these conditions and working modulo p^, there are (p'^ ~P^~^) such 
choices for each, and of the resulting (04, cg) pairs we noted above that p'^~^{p — 1) 
of them have p^\D. So, given a curve with p \ C4,ce, we have a probability of 
l//-i(p-l) that/|L>, which gives 1 /p'' for the probability that p''\\D. In symbols, 
we have that (for p > 5 and fc > 1) 

Prob p'^ II (c|-c^) ptc4,C6 =l/p''. 

Including the conditional probability for minimal models, we get 

(3) P™(p, fc) = (1 - l/p^°)-^(l - l/p)V/ (for P > 5 and fc > 1). 

Note that summing this over fc > 1 gives Kp(l — l/p)/p for the probability for an 
elliptic curve to have multiplicative reduction at p. 

What is the probability Pa(p, fc) that a curve given by a minimal model has 
additive reduction at p > 5 and p''\\D for some fc > 0? We shall temporarily ignore 
the factor of Kp = (1 — 1/p^")^^ from non-minimal models and include it at the 
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end. We must have that p\c4,ce, and thus get that k > 2. For k = 2,3,4, which 
correspond to Kodaira symbols II, III, and IV respectively, the computation is not 
too bad: we get that p^\\D exactly when p|c4 and p\\cq, so that the probability 
is (l/_p) • (1 — l/p)/p = (1 — l/p)/p^; for p^\\D we need p||c4 and p^jce and thus 
get (1 — \/p)/p ■ (1/p^) = (1 — ^/p)/p^ for the probability; and for we need 

J3^|c4 and p^||c6 and so get • (1 — l/p)/p^ = (1 ~ ^/p)lP^ for the probability. 

Note that the case A; = 5 cannot occur. Thus we have (for p > b) the formula 
Pa{p,k) = (1 - - for k = 2,3,4. 

More complications occur for fc > 6, where now we split into two cases depending 
upon whether additive reduction persists upon taking the quadratic twist by p. This 
occurs when p^|c4 and p'^\cq, and we denote by Pa{p, k) the probability that p''\\D 
in this subcase. Just as above, we get that Pa{p, k) = {1 — — \/p)lp^^^ 

for k = 8,9,10. These are respectively the cases of Kodaira symbols IV*, III*, 
and II*. For fc = 11 we have P"(p, fc) = 0, while for k > 12 our condition of 
minimality implies that we should take Pa{p, k) = 0. 

We denote by Pa{p, k) the probability that p^\D with either p^||c4 or p'^Ucg. First 
we consider curves for which p'^\D, and these have multiplicative reduction at p upon 
twisting. In particular, these curves have p^||c4 and p^||c6, and the probability of 
this is (1 - l/p)/p^ • (1 - l/p)/p^. Consider fc > 7. We then take ca/p^ and cq/p^ 
both modulo p^~^ ^ and get that p^~^\\{D /p^) with probability l/p^~^ in analogy 
with the above. So we get that P*{p,k) = (1 - 1/jj1°)-^(1 - Ijpflp^'^ for fc > 7. 
This corresponds to the case of I^.g- 

Finally, for p^\D (which is the case \%) wo get a probability of (1/p^) • (1/p^) 
for the chance that p^|c4 and p'^jcg, and (since there are p points mod p on the 
auxiliary curve {ci/p^)^ = {c^/p^Y (mod p)) a conditional probability of {p'^—p)/p'^ 
that p^\\D. So we get that P*(p, 6) = (1 - l/p^^)-^{l - l/p)/p^. 

We now impose our current notation on the previous paragraphs, and naturally 
let Pl{p, fc) = and P^{p, fc) = Pa{p, k) for fc < 5. Our final result is that 



(4) P^{p,k) = 



(5) Pl{p,k) 



L - l/piO)-^(l - 1/p)// fc = 2,3,4 

1 - l/pi°)-Hl - l/p)/^''"^ fc = 8,9,10 

'(1 - l/piO)-i(l - l/p)/p5 fc = 6 

(1 - l/piO)-i(l - 1/p)Vp''"^ fc>7 



with P^{p,k) and P*{p,k) equal to zero for other fc. We conclude by defining 
Po{p,k) to be zero for fc > and to be the probability (1 — l/p^°)~^(l — 1/p) 
that p I D for fc = 0. We can easily check that we really do have the required 
probability relation J2T=o [Pm{p, k) + P^{p, fc) + P*{p, fc) + Po(p, fc)] =1, as: the 
cases of multiplicative reduction give Kp{l — l/p)/p] the cases of Kodaira symbols 
II, III, and IV give Kp(l/p^ — 1/p^); the cases of Kodaira symbols IV*, III*, and II* 
give Kp{l/p'^ — l/p^°); the cases of I^ summed for fc > 1 give Kp{l — l/p)/p^; the 
case of Iq gives Kp{l — l/p)/p^] and the sum of these with Po{p,Q) = Kp{l — 1/p) 
does indeed give us 1. We could do a similar (more tedious) analysis for p = 2,3, 
but this would obscure our argument. 

Given a curve of discriminant D, we can now compute the expectation for its 
Tamagawa number. We consider primes p\D with p > 5, and compute the local 
Tamagawa number t{p). When E has multiplicative reduction at p and p'^H-D, then 
t{p) = k if — C6 is square mod p, and else t{p) = 1,2 depending upon whether fc is 
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odd or even. So the average of \/t{p) for this case is tm{k) = ^(1 + Vfc), \{\/^+Vk) 
for k odd/even respectively. 

When E has potentiaUy multiphcative reduction at p with p^\\D, for k odd 
we have t{p) = 4,2 depending on whether (cg/p^) ■ (A/p'^) is square mod p, and 
for k even we have t{p) = 4,2 depending on whether A/p'' is square mod p. In 
both cases the average of is i(^/2 + V4). In the case of Iq reduction where 

we have p^\\D, we have that t{p) = 1,2,4 corresponding to whether the cubic 

— (27c4/p^)a; — {54:Cq/p^) has 0, 1, 3 roots modulo p. So the average of \/t{p) is 

Vl{ip ~ l)ip + l)/3) + V2{p{p - l)/2) + Va{{p - l){p - 2)/6) _ 2 V2 1 

((p-l)(p+l)/3) + (p(p-l)/2) + ((p-l)(p- 2)/6) 3 2 3p- 

in this case. 

For the remaining cases, when p'^\\D or p^^'HD we have t{p) — 1, while when 
or p^ll-D we have t{p) — 2. Finally, when we have t{p) = 3, 1 depending 

on whether —6ce/p^ is square mod p, and similarly when p®||-D we have t{p) = 3, 1 
depending on whether — Bcg/p* is square mod p, so that the average of ■\/t{p) in 
both cases is i(l + V3). We get that e^ik) = 1, \/2, i(l + V3), ^(1 + V3), V2, 1 for 
,9, 10, while 

+ fcodd (1+^^^, 

i(V2 + \/fc), fceven ^ \i(V2 + \/4), fc>7 

with e"(fc) and e^(p. A:) equal to zero for other k. 

We define the expected square root of the Tamagawa number K{p) at p by 

oo 

(7) K{p) = ^ [em{k)P„,{p, k) + el{k)P:{p, k) + 6^(p, fc)P,*(p, fc) + P^{p, k)] 

k=0 

and the expected global^ Tamagawa number to be (3r — Y[p K{p). The convergence 
of this product follows from an analysis of the dominant k — 0, 1,2 terms of (O, 
which gives a behaviour of 1 + 0(l/p^). So we get that the Tamagawa product is a 
constant on average, which we do not bother to compute explicitly (we would need 
to consider p = 2, 3 more carefully to get a precise value). 

To compute the average value of aA{E) = Yip F{p) in HI we similarly assume'^ 
that each prime acts independently; we then compute the average value for each 
prime by calculating the distribution of F{p) when considering all the curves mod- 
ulo p (including those with singular reduction, and again making the slight ad- 
justment for non- minimal models). This gives some constant for the average a a 
of aA{E), which we do not compute explicitly. Note that Jlp-^l-P) converges if 
we assume the Sato- Tate conjecture j31) since in this case we have that is p on 
average. 

4. Relation between conductor and discriminant 

We now give heuristics for how often we expect the ratio between the absolute 
discriminant and the conductor to be large. The main heuristic we derive is: 



''Note that the Tamagawa number at infinity is 1 when E has negative discriminant and else 
is 2, the former occurring approximately \/3/(l + VS) ~ 63.4% of the time. 

^This argumentative technique can also be used to bolster our assumption that using Connell's 
conditions should be independent of other considerations. 
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Heuristic 4.1. The number B{X) of rational elliptic curves whose conductor is 
less than X satisfies B{X) ^ cX^^^ for an explicit constant c > 0. 

To derive this heuristic, we estimate the proportion of curves with a given ratio 
of (absolute) discriminant to conductor. Since the conductor is often the squarefree 
kernel of the discriminant, by way of explanation we first consider the behaviour 
of f{n) — n/sqfree(n). The probability that f{n) = 1 is given by the probability 
that n is squarefree, which is classically known to be 1/C(2) = 6/7r^. Given a 
prime power p™, to have f{n) = says that n = p^^^u where u is squarefree 
and coprime to p. The probability that p™+^||ri is (1 — l/p)/p™+^, and given this, 
the conditional probability that is squarefree is (6/7r^) • (1 — l/p^)~^. 

Extending this multiplicatively beyond prime powers, we get that 

P Kf / f r ^ 1 6 -pr 6 1-p. 1 
Prob n/sqfree(n) = q\ = — — — - = \\ -. 

p"ll« p\q 

In particular, the average of f{n)"' exists for 7 < 1; in our elliptic curve analogue, 
we will require such an average for 7 = 5/6. We note that it appears to be an 
interesting open question to prove an asymptotic for ^ n/sqfree(n). 

n<X 

4.1. Derivation of the heuristic. We keep the notation Z) = |A| and wish to 
compute the probability that D/N = q for a fixed positive integer q. For a prime 
power p" with p > 5, the probability that p'"\\{D/N) is given by: the probability 
that E has multiplicative reduction at p and p'"^^\\D, that is Pm{p, v + 1); plus the 
probability that E has additive reduction at p and p'"^'^\\D, that is Pa{p,v + 2); 
and the contribution from Pq(p,v), which is zero for v > and for w = is the 
probability that p does not divide D. So, writing v — Vp{q), we get that (with a 
similar modified formula for p = 2,3) 

(8) Proh[D/N = q] ='[[Ep{vp{q))=l[[Prn{p,l + v)+Paip,2 + v)+Po{p,v)]. 

p p 

It should be emphasised that this probability is with respect to (as in the previous 
section) the ordering of the curves by discriminant. We have 
00 00 

(9) E 1 ^ E E P^oh[D/N = g] ^ ^ a{qXf/' ■ PT0h[D/N = q] , 

E:Ne<X q=lE:D<qX q=l 

where a — a+ + a_ from the Brumer-McGuinncss heuristic 12.11 If this last sum 
converges, then we get Heuristic 14.11 

To show the last sum in ^ does indeed converge, we upper-bound the probability 
in ©. We have that Pm{p,v + \) < l/p"+^ and Pa(p, v + 2) < 2/p"+^, which implies 



q p 

p\q 



f{q) = Proh[D/N = g] = [] Ep{vp{q)) < ^ [] 

p 

We then estimate 

f: f ^ n n +£ (^) s n (1 + 

q=l q=l ^ p\q ' p ^ 1=1 ^ ' p ^ ^ 

and the last product is convergent upon comparison to <^(7/6)'^. Thus we shown 
that the last sum in © converges, so that Heuristic 14 . 1 1 follows . 
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We can note that Fouvry, Nair, and Tencnbaum 13_ have shown that the number 
of minimal models with D < X is at least cX^^^, and that the number of curves 
with D < X with Szpiro ratio > k is no more than CeX^/'*+' for every e > 0. 

4.2. Dependence of D/N and the Tamagawa product. We expect that D/N 
should be independent of the real period, but the Tamagawa product and D/N 
should be somewhat related.® We compute the expected square root of the Tama- 
gawa product when D/N = q. As with © and using the e defined in we find 
that this is given by 

. , T-r [emivi)Pmip,vi) + e:iv2)P^ip,V2) + eiip,V2)P^ip,V2) + Poip,v)] 
\^ [P„,{p,v^)+Paip,V2) + Poip,v)] 

where vi = v + 1 and V2 = v + 2 and v =^ Vp{q). 

4.3. The comparison of log A with logiV. We now want to compare log A 
with logA^, and explicate the replacement therein in Guess In order to bound 
the effect of curves with large D/N, we note that 

Proh[D/N>Y]^Y.f(l)^Y.l]ll^ 

q>Y q>Y ^ p\q ^ 

and use Rankin's trick, so that for any < a < 1 we have (using — 1 > alogp) 

Prob[l./A. > r] - ^n(l + J:^ + + • 

Y^rif 3/p \ /m-^ c/p \ e'=^'^i°s^ 



1 I ^ V 

a logp / Y 



for some constants c, c, by taking a = I /^/logY (this result is stronger than needed) . 

However, a more pedantic derivation of Guess lX^ does not simply allow replacing 
logA^ by log A, but requires analysis (assuming il.^c{E) to be independent of q) of 



anaA 



3456 C(10) 



J I n,,iE)- ^r,(g)(logA/g)3/«.Prob[i?/iV = gj 



The above estimate on the tail of the probability and a simple bound on ri{q) in 
terms of the divisor function shows that we can truncate the q-sum at Y with an 
error of 0{1/Y^^^), and choosing (say) Y = e^'°s^ gives us that \og{A/q) ~ log A 
(note that we restricted to A > Vx) . So the bracketed term becomes the desired 

^(logA)3/8r;(g) • Pioh[D/N g] - /3^(log A)^/®, 

q<Y 

upon noting that the g-part of the sum converges to Pr as K — > cxd. 



The size of the torsion subgroup should also be related to D/N, but in the next section we 
argue that curves with nontrivial torsion are sufficiently sparse so as to be ignored. 
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4.4. Counting curves with vanishing L-value. We now estimate the number 
of elliptic curves E with even parity and L{E, 1) = when ordered by conductor. 

Heuristic 4.2. Let R{X) be the number of elliptic curves E with even parity and 
conductor less than X and L{E, 1) = 0. Then R{X) - cX'^'^/'^^{\ogXf/^ for some 
explicit constant c > 0. 

From Guess 1^21 we get that the number of even parity curves with < A < qX 
and D/N = q and L{E^ 1) = is given by 

W{qX) ■ {il{q)/(3r) ■ Prob[l?/7V = q] , 

and we sum this over all q. As we argued above, the tail of the sum does not affect 
the asymptotic (and so we can take log A ^ logA^ in W), and again we get that 
the g-sum converges. This then gives the desired asymptotic for the number of 
even parity curves with conductor less than X and vanishing central L-value (upon 
arguing similarly for curves with negative discriminant). 

5. Torsion and isogenies 

We can also count curves that have a given torsion group or isogeny structure. 
For instance, an elliptic curve with a 2-torsion point can be written as an integral 
model in the form = + ax^ + bx where A = 166^ (a^ — 46); thus, by lattice-point 
counting, we estimate about ^/X curves with absolute discriminant less than X. 
The effect on the conductor can perhaps more easily be seen by using the Fricke 
parametrisation ca = {t + 16) (t + 64)r^ and cq = {t - 8){t + 64)2T^ of curves 
with a rational 2-isogeny, and then substituting t — p/q and V = T/q to get 
C4 = {p+l6q)(j,+6Aq)V'^ and cq = {p-8q){p+6Aq)'^V^ so that A = p(j)+6Aqfq'^V^ . 
The summation over the twisting parameter V just multiplies our estimate by a 
constant, while ABC-estimates imply that there should be no more than X^^^^'^ 
coprime pairs {p, q) with the squarefree kernel of pq{p + 64q) smaller than X in 
absolute value. So we get the heuristic that almost all curves have no 2-torsion, 
even under ordering by conductor. Indeed, the exceptional set is so sparse that we 
can ignore it in our calculations. A similar argument applies for other isogenies, and 
more generally for splitting of division polynomials. Also, the results of Duke |12| 
for exceptional primes are applicable here, albeit with a different ordering. 

6. Experiments 

We wish to provide some experimental data for the above heuristics. In par- 
ticular, the two large datasets of Brumer-McGuinness [3] and Stein- Watkins PU] 
for curves of prime conductor up to 10^ and 10^° show little drop in the observed 
average (analytic) rank. Brumer and McGuinness considered about 310700 curves 
with prime conductor less than 10^ and found an average rank of about 0.978, while 
Stein and Watkins extended this to over 1 1 million curves with prime conductor up 
to 10^° and found an average rank of about 0.964. Both datasets are expected to 
be nearly exhaustive^ amongst curves with prime conductor up to the given limit. 
These results led some to speculate that the average rank might be constant. To 
test this, we chose a selection of curves with prime conductor of size 10^''. It is 
non-trivial to get a good dataset, since we must account for congruence conditions 
on the elliptic curve coefficients and the variation of the size of the real period. 



is one reason to take prime conductor curves; we also have | A| = with few exceptions. 



12 



MARK WATKINS 



6.1. Average analytic rank for curves with prime conductor near lO^''. 

As in , we divided the (c4,C6) pairs into 288 congruence classes witli (c4,C6) — 
(c4 mod 576, cg mod 1728). Many of tliese classes force the prime 2 to divide the 
discriminant, and thus do not produce any curves of prime conductor. For each 
class (64, cg), we took the 10000 parameter selections 

(C4,C6) = (576(1000 + i) + 54,1728(100000 + + 56) for (ij) G [1..10] x [1..1000], 

and then of these 2880000 curves, took the 89913 models that had prime discrimi- 
nant (note that all the discriminants are positive). This gives us good distribution 
across congruence classes, and while the real period does not vary as much as 
possible, below we will attempt to understand how this affects the average rank. 

It then took a few months to compute the (suspected) analytic ranks for these 
curves. We got about 0.937 for the average rank. We then did a similar experiment 
for curves with negative discriminant given by 

(C4,C6) = (576(-883 + i) +54,1728(100000 + j) + cq) for (i,j) e [1..10] x [1..1000], 

took the subset of 89749 curves with prime conductor, and found the average rank 
to be about 0.869. This discrepancy between positive and negative discriminant is 
also in the Brumer-McGuinness and Stein- Watkins datasets, and indeed was noted 
in [3].^'' We do not average the results from positive and negative discriminant; the 
Brumer-McGuinness Conjecture 12.11 implies that the split is not 50-50. 

In any case, our results show a substantial drop in the average rank, which, at 
the very least, indicates that the average rank is not constant. The alternative 
statistic of frequency of positive rank for curves with even parity also showed a 
significant drop. For prime positive discriminant curves it was 44.1% for Brumer- 
McGuinness and 41.7% for Stein- Watkins, but only 36.0% for our dataset — for 
negative discriminant curves, these numbers are 37.7%, 36.4%, and 31.3%. 

6.2. Variation of real period. Our random sampling of curves with prime con- 
ductor of size 10^* must account for various properties of the curves if our results 
are to possess legitimacy. Above we speculated that the real period plays the most 
significant role, and so we wish to understand how our choice has affected it. 

To judge the effect that variation of the real period might have, we did some 
comparisons with the Stein- Watkins database. First consider curves of positive 
prime discriminant, and write E as = Ax^ + b2x'^ + 2640; + 65 and ei > 62 > 63 for 
the real roots of the cubic. We looked at curves with even parity and considered the 
frequency of positive rank as a function of the root quotient t = f^^Eff , noting that"'^'^ 
^ 2'/'>7v{t-ty ^<' ^ r^^^ curves we considered all had 0.617 < t < 0.629. 

However, similar to when we considered curves ordered by conductor, before 
counting curves with extra rank, we should first simply count curves. Figure ^ 
indicates the distribution of root quotient t for the curves of prime (positive) dis- 
criminant and even parity from the Stein- Watkins database (more than 2 million 
curves meet the criteria). The x-axis is divided up into bins of size 1/1000; there are 
more than 100 times as many curves with t < 0.001 as there with 0.500 < t < 0.501, 
with the most extremal dots not even appearing on the graph. 

^"^"An interesting phenomenon was the systematic influence of the discriminant sign on all 
aspects of the arithmetic of the curve." 

^^The calculation follows as in the previous sections; via calculus, we can compute that this 
function is maximised at t ^ .0388505246188 with a maximum just below 4.414499094. 
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Next we plot the frequency of L{E, 1) = as a function of the root quotient in 
FigureEl Since there are only about 1000 curves in some of our bins, we do not get 
such a nice graph. Note that the left-most and especially the rightmost dots are 
much below their nearest neighbors, the graph slopes down in general, and drops 
more at the end. We see no evidence that our results should be overly biased. In 
particular, the frequency of L{E, 1) = is 41.7% amongst all even parity curves of 
prime discriminant in the Stein- Watkins database, and is 42.8% for the 12324 such 
curves with 0.617 < t < 0.629. The function plotted (labelled on the right axis) in 



Figure His of Q.^A^/'^^ 



agm(l,yt) 



as a function of t, and note that this goes 



to as t — *■ 0, 1; there is nothing canonical about the choice of our t-parameter, 
and we chose it more for convenience than anything else. 

Similar computations can be made in the case of negative discriminant, which we 
briefly discuss for completeness (again restricting to curves with even parity where 
appropriate). Let r be the real root of the cubic polynomial Ax^ + 62^^ -I- 264a; + 6g, 
and Z > the imaginary part of the conjugate pair of nonreal roots. Letting 
f — r + 62/12 and c — f/Z we then have^^ 

7rV2 



(1 



9c2/4)i/i2agm(l, 



3c 



4^1+9c2/4y 

We renormalise via taking C — 1/2 + arctan(c)/7r, and graph the distribution of 
curves versus C in Figure |2| The symmetry ^'^ of the graph might indicate that the 
coordinate transform is reasonable. All our curves have 0.555 < C < 0.557. 
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Figure 3. A < 0: Distribution of curves as a function of C 



^^This is maximised at c Ki —33.58515148525, with the maximum a bit less than 8.82921518. 

"'^'^The blotches around 0.22-0.23 and 0.77-0.78 appear to come from the fact that curves with 
a4 small (in particular ±1) tend to have C in these ranges (for our discriminant range), and this 
causes instability in the counting function. 
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Next we plot the frequency of L{E, 1) = as a function of the root quotient in 
Figure^! Again we also graph the function rirc|A|^/^^ on the right axis. Here the 
drop-off is more pronounced than with the curves of positive discriminant. Note 
the floating dot around C = 1/2. Indeed the 100 closest curves with C < 1/2 all 
have positive rank; this breaks down when crossing the 1/2-barrier. This is not 
particularly a mystery; these curves have ae — and/or be = 1, and thus have an 
obvious rational point. Recall that C = 1/2 corresponds to c = = f. 



1 I 1 1 1 1 1 9 



0.8 - - 7.2 



0.6 - - 5.4 




I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

Figure 4. A < 0: Positive rank frequency as a function of C, 
and ilrclAI^/^^ as a function of C. 

We again see no evidence that our results should be biased. In particular, the 
frequency of L{E^ 1) = is 36.4% amongst all even parity curves of negative prime 
discriminant in the Stein- Watkins database, and is 37.0% for the 4695 such curves 
with 0.555 < C < 0.557. 

6.3. Other considerations. The idea that the "probability" that an even parity 
curve possesses positive rank should be proportional to V^rc is perhaps overly 
simplistic — in particular, it is not borne out too precisely by the Stein- Watkins 
dataset. We consider positive prime discriminant curves with even parity; for those 
with 0.64 < fire < 0.65 we have 78784 curves of which 45.9% have positive rank, 
while of the 9872 with 0.32 < ili-o < 0.325 we have 36.0% with positive rank, for 
a ratio of 1.28, which is not too close to v^- One consideration here is that we 
have placed a discriminant limit on our curves, and there are curves with larger 
discriminant and 0.32 < ili-c < 0.325 that we have not considered. This, however, 
is extra-particular to the idea that only the real period should be of import. 

One possibility is that curves with small discriminant and/or large real period 
have smaller probability of L{E, 1) = that our estimate of cy/Tl^ would suggest — 
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indeed, it might be argued (maybe due to arithmetic considerations, or perhaps ex- 
phcit formulae for the zeros of L- functions) that curves with such smaU discriminant 
cannot realise their nominal expected frequency of positive rank. Unfortunately, 
we cannot do much to quantify these musings, as the effect would likely be in a 
secondary term, making it difficult to detect experimentally. Note also that a rel- 
ative depression of rank for small discriminant curves would give a reason for the 
near-constant average rank observed by Brumer-McGuinness and Stein- Watkins. 

6.4. Mordell-Weil lattice distribution for rank 2 curves. We have other 
evidence that curves of small discriminant might not behave quite as expected. 
We undertook to compute generators for the Mordell-Weil group for all 2143079 
curves of (analytic) rank 2 of prime conductor less than 10^° in the Stein- Watkins 
database.^"* J. E. Cremona ran his mwrank programme |^ on all these curves, and 
it was successful in provably finding the Mordell-Weil group for 2114188 of these. 
For about 2500 curves, the search region was too big to find the 2-covering quartics 
via invariant methods, while around 8500 curves had a generator of large height 
that could not be found, and over 18000 had 2-Selmer rank greater than 2. We 
then used the FourDescent machinery of MAGMA [2] which reduced the number 
of problematic curves to 54. Of these, 19 have analytic III of 16.0 and we expect 
that either 3-descent or 8-descent will complete (assuming GRH to compute 
the class group) the Mordell-Weil group verification; for the 35 other curves, there 
is likely a generator of height more than 225 which we did not attempt to find."'^^ 

We then looked at the distribution of the Mordell-Weil lattices obtained from 
the induced inner product from the height pairing; since all of our curves have 
rank 2, we get 2-dimensional lattices. We are not so interested in the size of the 
obtained lattices, but more so in their shape. Via the use of lattice reduction (which 
reduces to continued fractions in this case), given any two generators we can find 
the point P of smallest positive height on the curve. By normalising P to be the 
unit vector, we then get a vector in the upper-half-plane corresponding to another 
generator Q. Via the standard reduction algorithm, we can translate Q so that it 
corresponds to a point in the fundamental domain for the action of SL2(Z). Finally, 
by replacing Q by —Q if necessary, we can ensure that this point is in the right half 
of the fundamental domain (in other words, we must choose an embedding for our 
Mordell-Weil lattice). In this manner, for each rank 2 curve we associate a unique 
point z = X + iy in the upper-half-plane with + > 1 and < a: < 1/2. 

With no other guidance, we might expect that the obtained distribution for the z 
is given by^^ the Haar measure {dxdy)/y'^ . We find, however, that this is not borne 

1/2 

out too well by experiment. In particular, we should expect that ~ 95.5% 
of the curves should have y > 1, while the experimental result is about 93.5%. 
Furthermore, we should expect that the proportion of curves with y > Y should 
die off like \/Y as y ^ oo; however, we get that 35.4% of the curves have y > 2, 



We also computed the Mordell-Weil group for curves with higher ranks but do not describe 
the obtained data here. 

bit more searching might resolve a few of the outstanding cases, but the extremal case of 
[0, 0, 1, -237882589, -1412186639384] appears to have a generator of height more than 600, and 
thus other methods will likely be needed to try to find it. T. A. Fisher has recently used 6-descent 
to find some of the missing points. 

^^Siegel 1271 similarly uses Haar measure to put a natural measure on n-dimensional lattices 
of determinant 1. 
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only 9.4% have y > 4, while 1.7% have y > 8 and 0.2% have y > 16. The vahdity of 
the vertical distribution data might be arguable based upon concerns regarding the 
discriminant cutoff of our dataset, but the horizontal distribution is also skewed. If 
we consider only curves with y > 1, then we should get uniform distribution in the 
x-aspect; however, Tabled shows that we do not have such uniformity. 

Table 1. Horiztonal distribution of rank 2 lattices with y > 1 



0.00 < 


X 


< 0.05 


9.0% 


0.25 < 


X 


< 0.30 


10.0% 


0.05 < 


X 


< 0.10 


9.6% 


0.30 < 


X 


< 0.35 


10.2% 


0.10 < 


X 


< 0.15 


9.8% 


0.35 < 


X 


< 0.40 


10.5% 


0.15 < 


X 


< 0.20 


9.9% 


0.40 < 


X 


< 0.45 


10.6% 


0.20 < 


X 


< 0.25 


10.0% 


0.45 < 


X 


< 0.50 


10.5% 



We cannot say whether these unexpected results from the experimental data are 
artifacts of choosing curves with small discriminant; it is just as probable that our 
Haar- measure hypothesis concerning the lattice distribution is simply incorrect. 

6.5. Symmetric power L-functions. Similar to questions about the vanishing 
of L{E,s), we can ask about the vanishing of the symmetric power L-functions 
L{Sym^^~^E, s). We refer the reader to ^ for more details about this, but mention 
that, due to conjectures of Deligne and more generally Bloch and Beilinson we 
expect that we should have a formula similar to that of Birch and Swinnerton-Dyer, 

stating that L{SyTc?^~^E, k){2'KN)^*'^) ' ^n^^^ should be rational with smaU de- 
nominator. Here, for k odd, 0+ is the real period and the imaginary period, with 
this reversed for k even. Ignoring the contribution from the conductor, and crudely 
estimating that Vlj^ « « l/A^^^^, an application of discretisation as before gives 
that the probability that L(Sym^'''^^i?, s) has even parity and L{SyTc?^^^ E , k) — 
is bounded above (cf. the ignoring of N) by c(log A)^/^- -y/l/A*^^^. Again following 
the analogy of above, we can then upper-bound the number of curves with conductor 
less than X with even-signed symmetric {2k — l)st power and L(Sym^'^~^i?, k) = 
by Cfc(e)X^/^^'' for every e > 0. It could be argued that we should order 

curves according to the conductor of the symmetric power L-function rather than 
that of the curve, but we do not think such concerns are that relevant to our im- 
precise discussion. In particular, the above estimate predicts that there are finitely 
many curves with extra vanishing when fc > 5. It should be said that this heuristic 
will likely mislead us about curves with complex multiplication, for which the sym- 
metric power L-function factors (it is imprimitive in the sense of the Selberg class), 
with each factor having a 50% chance of having odd parity. However, even ignoring 
CM curves, the data of (22j find a handful of curves for which the 9th, 11th and 
even the 13th symmetric powers appear (to 12 digits of precision) to have a central 
zero of order 2. We find this surprising, and casts some doubt about the validity of 
our methodology of modelling of vanishings. 

6.6. Quadratic twists of higher symmetric powers. The techniques we used 
earlier in this paper have also been used to model vanishings in quadratic twist 
families, and we can extend the analyses to symmetric powers. 
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6.6.1. Non-CM curves. We fix a non-CM curve E and let Ed be its dth qua- 
dratic twist, taking d to be a fundamental discriminant. From an analogue of 
the Birch-Swinnerton-Dyer conjecture we expect to get a small-denominator ra- 
tional from the quotient^"^ L{SyTn^ Ed,2){2T:NE) /^\m{Ed)^VL^c{Ed)- We have that 
^im{EdY^rc{Ed) ~ ^u-aiE)/ d?^"^ ' ^rc{E) / ^ and so we expect the number of 
fundamental discriminants \d\ < D such that L(S'ym^Ed, s) has even parity with 
L{Sym^Ed, 2) = to be given crudely (up to log-factors) by J2d<D c/V(P- So we 
expect about (logZ?)^ quadratic twists with double zeros for the 3rd symmetric 
power; generalising predicts finitely many extra vanishings for higher (odd) powers. 



Table 2. Fundamental d with ordL(Sym^i?(i, s) > 2 



11a 


-40 -52 -563 -824 -1007 -1239 -1460 -1668 -1799 -2207 
-2595 -2724 -2980 -3108 -3592 -4164 -4215 -4351 -4399 
12 69 152 181 232 273 364 401 412 421 444 476 488 652 669 696 
933 1101 1149 1401 1576 1596 1676 1884 1928 2348 2445 2616 
2632 3228 3293 3404 3720^ 3793 4060 4093 4161 4481 4665 4953 


14a 


-31 -52 -67 -87 -91 -111 -203 -223 -255 -264 -271 -311 -327 
-367 -535 -552 -651 -759 -804 -831 -851 -852 -920 -1099 -1263 
-1267 -1335 -1524 -1547 -1567 -1623 -1679 -1707* -2047 -2235 
-2280 -2407 -2443 -2563 -2824 -2831 -3127 -3135 -3523 -4119 
-4179 -4191 -4323 

137 141* 229 233 281 345 469 473 492 497* 697 901 1065 1068 1353 1457 
1481 1513 1537 1793 1873 1905 2024 2093 2193 2265 2321 2589 2657 2668 
2732 2921 2981 2993 3437 3473 3529 4001 4124 4389 4488 4661 4817 


15a 


-11 -51 -71 -164 -219 -232 -292 -295 -323* -340 -356 -399 -519 
-580 -583 -584 -671 -763 -804 -851 -879 -943 -1012 -1060 -1151 
-1199 -1284 -1288 -1363 -1551 -1615 -1723 -1732 -2279 -2291 
-2379 -2395 -2407 -2571 -2632 -2635 -2756 -3396 -3588* -3832 
17 21 61 77 136 156 181 229 349 444 481 501 545 589 649 781 876 905 
924 949 1009 1144 1249 1441 1501 1580 1621 1804 1861 1921 2041 2089 
2109 2329 2581 2829 2840 2933 3001 3916 



We took the curves lla:[0, —1, 1, 0, 0] and 14a:[l, 0, 1, —1, 0], and computed either 
LlSyui'^ Ed, 2) or L'(Sym^Ed, 2) for all fundamental discriminant d with \d\ < 5000. 
We did the same for 15a: [1, 1, 1, 0, 0] for |d| < 4000. We then looked at the number 
of vanishings (to 9 digits of precision). For 11a we found 58 double zeros and one 
triple zero (indicated by a star in Table while for 14a we found 88 double zeros 
and three triple zeros, and 15a yielded 83 double zeros and two triple zeros. 

6.6.2. CM curves. Next we consider CM curves, for which we can compute signifi- 
cantly more data, but the modelling of vanishings is slightly different. Let E he a 
rational elliptic curve with CM, and ip its Hecke character. We shall take ip to be 
"twist-minimal" — this is not the same as the "canonical" character of Rohrlich , 
but rather we just take i? to be a minimal (quadratic) twist. Indeed, we shall only 

contribution from the conductor actually comes from non-integral Tamagawa numbers 
from the Bloch-Kato exponential map, and in the case of quadratic twists, the twisting parameter d 
should not appear in the final expression. 
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consider 11 different choices of E, given (up to isogeny class) by 27a, 32a, 36a, 
49a, 121a, 256a, 256b, 361a, 1849a, 4489a, and 26569a, noting that 27a/36a and 
32a/256b are respectively cubic and quartic twist-pairs. In our tables, these can 
appear in a briefer format, such as 67^ for 4489a. 

We normalise the Hecke i-function L{ip, s) to have s = 1 be the center of 
the critical strip. For d a fundamental discriminant, we let ipd be the Hecke 
Grossencharacter -0 twisted by the quadratic Dirichlet character of conductor d. 
Finally, note that the symmetric powers L{Sym''^p, s) are just L{'4}^,s), where we 
must take ip^ to be the primitive underlying Grossencharacter. 

We then expect L(^'^, 2)(27r)/riim(i5)'^ to be rational with small denominator. 
We can then use discretisation as before to count the expected number of funda- 
mental discriminants \d\ < D for which the L-function L^Tp'^, s) has even parity but 
vanishes at the central point — since we have flimiEd)"^ ~ flimiE) / d'^^'^ , we expect 
the number of discriminants d that yield even parity and L{ip^,2) = is crudely 
given by J2d<D l/v(^^, so we should get about D^/'^ such discriminants up to D. 

For higher symmetric powers, we expect that L('0^'^~^, A:)(27r)''~^/ri_|_(i5)^'^~^ is 
rational with small denominator, and thus get that there should be finitely many 
quadratic twists of even parity with vanishing central value. 

We took the above eleven CM curves and took their (fundamental) quadratic 
twists up to 10^. We must be careful to exclude twists that are isogenous to other 
twists. In particular, we need to define a primitive discriminant for a curve with 
CM by an order of the field K — this is a fundamental discriminant d such that 
disc(if) does not divide d, expect for if = Q(i) when d > is additionally primitive 
when 8\\d. Note also that 27a and 36a have the same symmetric cube L-function. 



Table 3. Counts of double order zeros for primitive twists 





27a 


32a 


36a 


49a 


121a 


256a 


256b 


361a 


1849a 


4489a 


26569a 


3rd 


59 


32 




67 


78 


32 


21 


45 


28 


31 


1 


5th 


3 


1 


5 


2 


1 


2 


2 














7th 








2 





1 





















TableOlists the following results for counts of central double zeros (to 32 digits) 
for the L-functions of the 3rd, 5th, and 7th symmetric powers.^* Tables 0] and |S1 list 
the primitive discriminants that yield the double zeros. The notable signedness can 
be explained via the sign of the functional equation. We are unable to explain the 
paucity of double zeros for twists of 26569a; Liu and Xu have the latest results PU] 
on the vanishing of such L- functions, but their bounds are far from the observed 
data. Similarly, the last-listed double zero for 4489a at 67260 seems quite small. 

There appear to be implications vis-a-vis higher vanishings in some cases; for 
instance, except for 27a, in the thirteen cases that L(V'^,s) has a double zero at 
s = 3 then L{'4>d,s) also has a double zero at s = 1. Similarly, the 7th symmetric 
power for the 27365th twist of 121a has a double zero, as does the 3rd symmetric 
power, while the L-function of the twist itself has a triple zero. Also, the 22909th 
twist of 36a has double zeros for its first, third, and fifth powers (note that 36a 
does not appear in Table 01 as the data are identical to that for 27a). 



"'^^We found no even twists with L(^^, 5) = 0, and no triple zeros appeared in the data. 
The local signs at p = 2, 3 involve wild ramification are more complicated (see EUEIE] 
for a theoretical description), and thus there is no complete correlation in some cases. 
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Table 4. Primitive d with ord s) = 2 



27a 



172 524 1292 1564 1793 3016 4169 4648 6508 9149 9452 9560 10636 
11137 12040 13784 14284 15713 17485 17884 22841 22909 22936 25729 
27065 27628 29165 30392 34220 35749 38636 40108 41756 44221 47260 
51512 54385 57548 58933 58936 58984 59836 59996 62353 64268 70253 
74305 77320 77672 78572 84616 86609 86812 87013 92057 95861 96556 
97237 99817 



32a 



-395 -5115 -17803 -25987 -58123 
-60347 -73635 -79779 -84651 -99619 
257 1217 2201 2465 14585 26265 45201 82945 
4632 5336 5720 7480 9560 30328 30360 
31832 38936 45848 69784 71832 83512 92312 



49a 



-79 -311 -319 -516 -856 -1007 -1039 -1243 -1391 -1507 -1795 
-2024 -2392 -2756 -2923 -3527 -3624 -4087 -4371 -4583 -4727 
-5431 -5524 -5627 -6740 -7167 -7871 -8095 -8283 -10391 -10628 
-13407 -13656 -13780 -16980 -18091 -22499 -27579 -28596 -30083 
-30616 -32303 -32615 -36311 -36399 -38643 -39127 -40127 -42324 
-52863 -64031 -64399 -66091 -66776 -66967 -69647 -70376 -71455 
-72663 -73487 -73559 -77039 -84383 -90667 -91171 -98655 -98927 



11^ 



12 140 632 1160 1208 1308 1704 1884 2072 2136 2380 2693 2716 3045 
4120 4121 5052 5528 5673 5820 6572 7532 11053 11208 12277 12568 
12949 13884 14844 15465 16136 18588 18885 19020 19884 24060 25788 
27365 27597 28265 28668 29109 29573 32808 32828 35261 36552 37164 
38121 38297 44232 44873 49512 49765 50945 52392 54732 55708 56076 
56721 58460 59340 65564 66072 66833 71688 72968 79557 80040 80184 
83388 84504 84620 84945 86997 87576 92460 95241 



256a 



401 497 2513 3036 3813 6933 6941 9596 9932 11436 14721 17133 17309 
18469 21345 21749 26381 26933 28993 29973 30461 33740 51469 53084 

62556 63980 67721 69513 73868 76241 81164 87697 



256b 



73 345 3521 5133 6693 7293 21752 25437 27113 34657 38485 41656 
42433 44088 46045 75581 79205 83480 89737 93624 96193 



19^ 



44 60 1429 1793 3297 3340 3532 3837 3880 4109 5228 5628 7761 8808 
9080 9388 12280 12313 12545 13373 13516 13897 19164 22204 23241 
25036 25653 41205 41480 42665 43429 44121 44285 44508 45660 48828 

50584 52989 64037 74585 75324 76921 81885 85036 96220 



43^ 



88 152 440 2044 4268 5852 6376 7880 8908 9880 14252 
15681 17864 20085 20353 28492 29477 45368 55948 56172 
57409 60177 68136 79916 84524 85580 86853 96216 



67^ 



17 57 869 1612 1628 3260 6380 6385 7469 8328 11017 13772 
14152 14268 14552 15901 22513 24605 24664 27992 29676 33541 
33789 36344 36588 38028 40280 43041 49884 62353 67260 



163^ 



30720 





Table 5. Primitive d with ord L{ipl'' 

s=k 




2 for some A; > 3 


27a 


5th: -13091 4040 18044 


49a 


5th: 437 19317 


32a 


5th: 1704 


121a 


5th: -183 7th: 27365 


36a 


5th: -856 -2104 -31592 -88580 22909 


256a 


5th: -79 -21252 


36a 


7th: -95 2488 


256b 


5th: -511 89320 
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6.6.3. Comparison between the CM and non-CM cases. For the twist computations 
for the symmetric powers, we can go much further (about 20 times as far) in the 
CM case because the conductors do not grow as rapidly.^° For the 3rd symmetric 
power, the crude prediction is that we should have (asymptotically) many more 
extra vanishings for twists in the CM case than in the non-CM case, but this is 
not borne out by the data. Additionally, we have no triple zeros in the CM case 
(where the dataset is almost 100 times as large), while we already have six for the 
non-CM curves. This is directly antithetical to our suspicion that there should be 
more extra vanishings in the CM case. As before, this might cast some doubt on 
our methodology of modelling of vanishings. 
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